Dynamical Casimir Effect for TE and TM Modes in a 
Resonant Cavity Bisected by a Plasma Sheet 
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Parametric photon creation via the dynamical Casimir effect (DCE) is evaluated numerically, in a 
three-dimensional rectangular resonant cavity bisected by a semiconductor diaphragm (SD), which is 
irradiated by a pulsed laser with frequency of GHz order. The aim of this paper is to determine some 
of the optimum conditions required to detect DCE photons relevant to a novel experimental detection 
system. We expand upon the thin plasma sheet model [Crocce et al., Phys. Rev. A 70 033811 
(2004)] to estimate the number of photons for both TE and TM modes at any given SD position. 
Numerical calculations are performed considering up to 51 inter-mode couplings by varying the SD 
location, driving period and laser power without any perturbations. It is found that the number of 
photons created for TE modes strongly depends on SD position, where the strongest enhancement 
occurs at the midpoint (not near the cavity wall); while TM modes have weak dependence on SD 
position. Another important finding is the fact that significant photon production for TMm modes 
still takes place at the midpoint even for a low laser power of 0.01 /ij/pulse, although the number 
of TEiii photons decreases almost proportionately with laser power. We also find a relatively wide 
tuning range for both TE and TM modes that is correlated with the frequency variation of the 
instantaneous mode functions caused by the interaction between the cavity photons and conduction 
electrons in the SD excited by a pulsed laser. 

PACS numbers: 42.50.Dv; 42.50.Lc; 42. 60. Da; 42.65.Yj 
Keywords: Cavity QED; Dynamical Casimir efEect 



I. INTRODUCTION 



Motion induced radiation, or the dynamical Casimir ef- 
fect (DCE) as it is now more commonly known, was first 
discussed by Moore [l| in 1970 who showed that pho- 
tons would be created in a Fabry-Perot cavity if one of 
the ends of the cavity wall moved with periodic motion. 
The mechanical oscillation amplitude of the wall is small: 
v/c 1, where v is the wall velocity and c is the speed 
of light and in this limit the number of photons produced 
during a given number of parametric oscillations is pro- 
portional to sinh^(2wt d/c), which was first discussed in 
Ei (for a nice review see, e.g. Q)- A realistic cavity 
has a finite quality factor, Q, and any exponential growth 
eventually saturates proportionally to smh^ {2Q v / c) , for 
times greater than the characteristic timescale r = Q/uj, 
e.g., see [1]. The emission of radiation from accelerated 
charges is a well known classical effect in electrodynam- 
ics. However, what is unusual with the DCE is the pre- 
diction that radiation can also be emitted from neutral 
moving objects. 

Aside from mechanical harmonic oscillations of a (cav- 
ity) wall [1, Q , it is also possible to induce temporal vari- 
ations of the dielectric function , also see [l3| . This 
leads to an effective wall motion by, for example, vary- 
ing the optical path length of the cavity [1, [ifl [ll] . A 
related approach to effective wall motion is to irradiate a 
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semiconductor sheet with a pulsed laser, which leads to a 
plasma mirror with time dependent surface conductivity, 
e.g., see [13]. (Note that pulsed lasers also have applica- 
tions in controlling the static Casimir force [isl.) This 
has been modeled by using thin dielectric slabs [H, HgI - 
[Tsj , or by using thin conducting slabs [l^ . Indeed exper- 
iments are already being built to detect DCE photons, 
based on the idea of usin g a semiconductor sheet irradi- 
ated by a pulsed laser [Iflim in centimeter sized cavities 
(for ideas relating to microcavities see e.g., [H, [2^.) In 
order to verify the DCE experimentally, it is important to 
consider some of the optimum conditions, which guides 
one to design and set up a detection system. 

In this paper we extend the model proposed by Crocce 
et al. [l^ to include both TE and TM modes for a three- 
dimensional resonant cavity (L^ x Ly x L^) bisected by a 
semiconductor diaphragm (SD). The SD is irradiated by 
a pulsed laser, which provides a periodic change between 
semiconducting and metallic states. Our primary con- 
cern is the dependence of the number of created photons 
upon the location of the SD, which can be performed only 
by numerical calculations, because perturbative treat- 
ments can no longer be applied to a general SD position 
away from the cavity wall at low temperatures. From 
an experimental point of view, it is desirable to divide 
the cavity into two parts, (i) a photon-creation/detection 
chamber and (ii) a pulsed-laser irradiation space to avoid 
the background caused by laser irradiation. Another im- 
portant factor to verify the DCE would be to keep the 
cavity at a low temperature of ^ 100 mK to suppress the 
number of thermal blackbody photons to less than unity. 
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Thus, numerical calculations are also performed to eval- 
uate the number of created photons under the conditions 
of high (50/xJ/pulse) and sufficiently low (0.01/iJ/pulse) 
laser powers. The results obtained are discussed in con- 
nection with our proposed DCE detection system (see 
Sec. irv|) using highly-excited Rydberg atom beams (Rb 
or K), which was already successfully applied to explore 
the dark matter axion 1241. 



II. FORMULATION 

One of the primary concerns of this work is the SD- 
position dependent enhancement of DCE photon cre- 
ation. Thus, we have extended the model of Crocce et 
al. [li], see Appendix El to both TE and TM modes for 
any SD position, rj ~ d/L^, where d is the distance of 
the SD from the cavity wall and is the longitudinal 
length of the cavity. Their work is closely related 
to the ID plasma sheet model of Barton & Calogeracos 
[2^. In Appendix Elwe give the 3D generalization, where 
Maxwell's equations are separated using scalar Hertz po- 
tentials |26l-l28l|. leading to two Klein-Gordon like equa- 



tions, given below. 

For TE modes a thin plasma sheet irradiated by a 
pulsed laser satisfies the following jump conditions (see 
Appendix |A|) at SD location z ~ d: 

disc*(r^,rf,t) 0, 

d\sc^'{v^,d,t) ^ ^^^>{r^,d,t) , (1) 

where e is the electron charge, ns{t) is the surface num- 
ber density of electrons excited to the conduction band, 
and m* is the effective mass of the conduction electrons 
(rj_ corresponds to transverse directions). Note that 
disc'^{d) = *(d+) - '^{d-). It is easy to show that the 
above jump conditions can be derived from the following 
self-adjoint differential equation 

Vi*(r, t) + *"(r, t) - *(r, t) = ^^^^S{z - d)*(r, t) , 

TO* 

(2) 

where 4'" and ^' represent derivatives with respect to z 
and t, respectively. The orthonormalized mode functions 
for a thin SD plasma sheet located at z = rf, then read 



A. 



(TE) 



B 



(TE) 



2 SUl{k,n,z) Wkj_(l"_L) , 

sm{kmSLz ~ z)) Vk^{r±) , 



< z < d 
d < z < Lz 



(3) 



where the TE transverse mode functions, iik^, are 




with V^Wkj^ = — k^Wk^. The case for TM modes is a 
little more delicate, but a short calculation leads to (see 
Appendix the jump conditions: 

disc<I)(r^,d,t) = f!j41$'(r^,d,t), 
disc^'{r±, d,t) = . (5) 
I 



These jump conditions can also be derived from the fol- 
lowing wave equation, 

Vi$(r,i)+$"(r,i)-$(r,<) = ^^S' {z - dMr,t) ; 

Klj^ TTl 

(6) 

however, this representation is not, strictly speaking, 
mathematically correct, e.g., sec [2^. The orthonormal- 
ized mode functions for TM modes are 



\ T^'^os{k,n^{Lz-z))u]^^{r±), d < z < 



I 

where the TM transverse mode functions are with V^Wkj^ = — k^^Wk^. Finally, by imposing continu- 

ity in conjunction with the jump conditions, we find the 

f X . ( TTm^x\ [2 . / TTmyy\ 



3 



1 e^n,{t) 



ki 



m* 



following eigenvalue relation: 
sin(fc,^^L^) 

(9) 

where the =F signs refer to TE and TM modes respec- 
tively (for TE drop the factor l/k"^). In the following 
we shall work with the dimensionless variable V{t)Lz ^ 
{e'^ns{t)/m*)Lz and thus, for TM modes the potential 
scales as V{t)Lz/(k.']_Ll). The time variation of km^{t) 
leads to a frequency variation in the angular frequency, 
'^m{t), see next section. 



of the SD then leads to ^ |3l|, ^ 

OO 



J2 M„eM^iP^ 



where the terms uJm{t) and 



\ dt 



is) 



(15) 



,*n , (16) 



III. QUANTIZATION 

In the plasma sheet model there is no time dependence 
in the bulk dielectric permittivity and so there are no is- 
sues with quantization in time-dependent dielectrics, e.g., 
see discussion in [l^ . Thus, quantization can be made 
straightforwardly with the quantum field operator expan- 
sion of the Hertz scalars [301 : 



(10) 



for TE modes, while expressions for TM modes are made 
by identifying ip ^ (j). The normalization is set to unity 
for reasons explained below. To evaluate the number 
of photons produced in a given mode we shall use the 
Bogolubov method by first definin g, f or t > (during 
irradiation) an instantaneous basis j3ll |: 



(11) 



with a similar expression for the TM component in terms 
of For TE modes with t < Q (before irradiation) we 
have 



V2< 



sinf ) Wk^(r±) (12) 



(for TM modes replace sin — > cos and — >■ ) where 
the stationary angular frequency is 



Ln 



(13) 

where c is the speed of light in vacuum. During irradia- 
tion {t > 0) the angular frequency becomes time depen- 
dent: 



correspond to squeezing and acceleration terms respec- 
tively (3^ . It is important to note that because of sym- 
metry, the relation for A/mn in the transverse (x, y) di- 
rections reduces to simple Kronecker deltas and hence, 
the mode functions reduce to an effective ID problem 
1331: the vector indices also become one-dimensional: 



m{mx,my,mz) 



This is also the reason 



why the normalisation is set to unity. As opposed to a 
vibrating cavity wall, simple sinusoidal expressions for 
Af,„„ do not exist, although for the plasma sheet model 
complicated analytic expressions can be found for a given 

■n- 

By defining auxiliary functions the problem reduces 
to a system of coupled first-order differential equations 
1331 1341 1 , also see M. These are solved by truncating the 
infinite sums with a cut-off in the mode sums such that 
^ (-max does not change the result in Eq. ((T9| below, 
as we have verified below. ^ The Bogolubov coefficients 
are defined as [33 



(v; 



rn 1 'r n 



/3„ 



-(v-rjcr) (17) 



where the invariant scalar product is (0, tj}) ~ 
—i /cavity '^^^('?^^* ~ Then, in terms of the origi- 

nal functions in Eqs. (|lip and (HI]), for out and in modes 
respectively, we find 



^™ p(n) 



1 



2a;,, 



in) 



(18) 

where amn is obtained by complex conjugation. The 
number of photons in a given mode (assuming an initial 
vacuum state) is given by^ 



(19) 



~T7 



(14) 



where {t) is the eigenvalue given in Eq. © . Substi- 
tution into the wave equation, ([2]) or ([B]), on either side 



^ The number of created photons saturates at a cut-ofT of imax 

50 in the mode sums. 
^ The number of photons created requires regularization such as 

by introducing an explicit frequency cut-off |32|| . 
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FIG. 1: (Color online) Example of pulse train (duration to), 
the laser pulse (period T) and asymmetric surface charge den- 
sity, ns{t), which is assumed as a flat-top plateau sandwiched 
by two asymmetric half Gaussians. 

where the total number of photons (if needed) is given 

by 



(20) 



An independent check can be made by confirming the 
unitarity condition which is satisfied by the Bogolubov 
coefficients: 



^(|a™„|'-|/3™„n = l, (21) 



where a„m and /3„m are defined above. 



IV. EXPERIMENTAL REMARKS 

Our primary concern is focused on cavity photons for 
TEiii and TMm modes, which have the lowest fre- 
quency where both TE and TM modes coincide. We shall 
consider an idealized lossless cavity with a rectangular 
shape {Lx = Lj, = Lr/2 = 0.05 m), which is bisected by 
a thin plasma sheet (SD) placed at z = d from one of the 
cavity walls. TEm and TMm modes have an angular 
frequency of wm = 3.05 x lO^^s^^ implying that we need 
a laser frequency of around ft = 2ajiii = 6.1 x 10^°[s~^] 
corresponding to a period of T = 111.1 [ps] (see Fig. [T]). 

In the present study, we assumed an infinitesimal thick- 
ness for the SD with a time-varying conductivity at a 
temperature below 1 K, where all the conduction elec- 
trons (those that are excited from the valence band and 
interact with the cavity photons) are expressed by the 
surface charge density (areal density) , Ug (t) (see Fig. [T]) . 
This is estimated assuming a laser wavelength of 860 nm 
and a penetration depth 5 /im for GaAs and is 

compatible with the diffusion length of the conduction 
electrons during laser-pulsed irradiation. In practice the 



thickness of the SD is about 0.5 mm and the penetration 
depth of the wm cavity photons in GaAs (carrier den- 
sity ^ lO^'^m^'^) when laser irradiated is estimated to be 
~ 50 ^m for a laser power of lOO/iJ/pulse. As discussed 
later, the plasma frequency ujp = yje'^ns/m*ds even for 
a low laser power of 0.01 /iJ/pulsc is much larger than 
^iii (by more than 10 times) and hence almost all the 
cavity photons are reflected by the plasma sheet. Thus, 
reflection takes place within a depth much less than the 
penetration depth in the SD when laser-irradiated. This 
is certainly small enough as compared with the scale of 
a microwave cavity of length ~ 100 mm. Taking into ac- 
count all these points, the assumption of an inflnitesimal 
thickness for the SD should be a good approximation to 
the real situation.'^ 

Another important factor is the relaxation time ((T2) to 
go from metallic to semiconducting states via recombina- 
tion immediately after pulse irradiation, where in general 
risit) has an asymmetric profile with a tail (see Fig. [T|). 
For unitary evolution the need for an asymmetric pro- 
file has been stressed by Dodonov & Dodonov[l6l fT7|: 
however, they also pointed out [l^ significant losses of 
photons in a semiconductor slab due to an asymmetric 
ns{t) profile with a long tail and thus, the need for a 
short relaxation time of less than about 20 ps to de- 
tect DCE photons with a frequency of 2.5 GHz. Re- 
cently, Agnesi et al. [2l| reported that the relaxation 
time could be reduced to ~ 20 ps for GaAs irradiated 
by neutrons. Alternatively, the relaxation time can be 
shortened by making a deep (trapping) level in the band 
gap of a semiconductor, which is introduced by high en- 
ergy electron irradiation and doping with Au etc. Indeed, 
it was reported that GaAs irradiated with 200 MeV Au"*" 
ions gives a relaxation time shorter than 4 ps at 7K [36} . 
Thus it is possible to make a profile without a long tail for 
V{t)Lz = e'^ns{t)Lz/Tn* which oscillates between and 
V-maxLz ranging typically from 1 to 10, 000 in our simu- 
lations. For the ns{t) profile we have assumed a flat-top 
plateau sandwiched by two asymmetric half Gaussians, 
as depicted in Fig. [1] with an off-set of = 35 ps, 
plateau oi tc — 7 ps and standard deviations of (Ti = 4 
and (72 = 8 ps, respectively. 

Much of the motivation for this work is based on a 
novel DCE detection system, using highly-excited Ryd- 
berg atom (Rb or K) beams working as a single photon 
detector, which has already been utilized to explore the 



^ Dodonov and Dodonov developed detailed analyses on the dis- 
sipation of cavity photons in a SD in a series of papers Il6l4l8ll . 
They discussed the temperature rise of a cavity caused by laser 
irradiation, which depends on the laser power (in their case as- 
sumed of order mj/pulse), thermal conductivity, and size of the 
SD. However, in our proposed experiment, the laser power will be 
constrained to a range of fij/ pulse with the duration of the pulse 
train set to the order of 100 ns (~ 1000 pulses) which is com- 
patible with the SD cooling time and thus, the temperature 
rise is likely to be very small. 
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FIG. 2: (Color online) Frequency variation for various /c„ for TE and TM modes with two different laser powers, in dimensionless 
units Vmaxiz = 5000 and Knaxiz = 1. Dashed lines denote baselines corresponding to stationary values. In Fig. [2] (c) the 
frequency variation is magnified by a factor of 10. 



dark matter axion [24|. The cavity is cooled down to 
100 niK by a dilute refrigerator using '^He combined 
with a pulse tube to suppress the number of thermal 
blackbody photons to be less than 0.5. Note that as well 
as DCE photons, thermal blackbody photons will also 
be enhanced parametrically [37j . so there is no a priori 
reason to do experiments at low temperatures. Our in- 
terest; however, is to detect the DCE as a pure vacuum 
fluctuation effect. At the same time, a low laser power 
is inevitable in order to keep the cavity at a low tem- 
perature and to achieve a high quality factor Q of better 
than 10^, which guarantees a lifetime of longer than a few 
milliseconds for the resonant microwave photons created. 

A high-sensitivity photon detector is required to ob- 
serve actual DCE photons and a highly excited Rydberg- 
atom beam is one of the most promising techniques, by 
working as a single photon detector in the microwave 
regime. Our concern in this work is to find the optimum 
conditions for parametric enhancement of DCE photon 
creation by numerical analysis and thus, we evaluate the 
number of photons created dependent on SD position, 
laser power, and the range of driving laser frequency giv- 
ing DCE photon creation. 

In general we must find solutions to fc„ (t) numerically 
in Eq. ©; however, in the limit 7y ^ 1 (or for low laser 
power) it is possible to find a linearized solution for TE 
modes only. This is similar to the approach of Crocce et 
al. [l^ who for ri = 1/2 with e'^ns{t)/m* ranging from 
10^° — 10^^ m~^, which is valid for doped semiconduc- 
tors at room temperature, also found a linearized solu- 
tion. Under practical experimental conditions; however, 
as mentioned, we require a low temperature of around 100 
mK to suppress thermal background radiation and to get 
a high Q value. If a laser power of 100 /xJ/pulse is also 
assumed, then the actual range is estimated to be from 
(non laser-irradiated) to around 10^ m~^ and hence, 
a perturbative analysis is not possible (unless 77 <C 1 for 
TE modes) . Indeed, this was part of our motivation for 



carrying out a full numerical and non-perturbative anal- 
ysis. 

It is also important to estimate numerically the fre- 
quency shift, which was recently discussed in [isl. [TqI. Issf . 
For the mechanical vibrations of a cavity wall the aver- 
age frequency (I)„ = k„c is close to that of the stationary 
angular frequency, wJJ, see Eq. (jl3p . However, for an SD 
irradiated by a pulsed laser this is not the case, because 
the time- variation of fc„(t) shifts from the baseline fre- 
quency asymmetrically (dashed lines in Fig. [2]) with a 
positive shift for TEii„ and a negative shift for TMii„ 
modes, see Fig. [5] for the 77 = 1/2 case. Hereafter, we 
denote uin — wJJ the "frequency shift," which should be 
taken into account when tuning the driving period of the 
pulsed laser for parametric photon creation. 



V. RESULTS & DISCUSSION 

The time variation of the eigenvalue kn{t) for both 
TEii„ and TMii„ modes is shown in Fig. ^ which re- 
flects the interaction between the cavity photons and the 
conduction electrons excited by the pulsed laser. Figs. 
[2] (a) and (b) show fc„(n = 1,2,3,4,5) for TE and TM 
modes, respectively, assuming rj = 1/2 and a Vma^Lz 
value of 5,000 (50^J/pulse). Here, we should note that 
the plasma frequency Wp ~ 3 x 10"^'^ s~^ is much larger 
than the frequency of the cavity photons of interestand 
thus, a pulsed laser will act as a high-frequency switch- 
ing (transmission/reflection) actuator. Laser irradiation 
leads to a positive change in A:„ for TEii„, while nega- 
tive changes occur for TMii„ modes from each baseline, 
which corresponds to the stationary fc^ value for the cav- 
ity photons. The width and variation of fc„(t) are almost 
independent of n for both TE and TM modes, but they 
gradually get smaller with increasing n. However, for TM 
modes the variation is significantly larger than that for 
TE modes, which is due to the strong coupling of the TM 
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Time Duration (ps) Time Duration (ps) 

FIG. 3: (Color online) Niii{t) with ri — 1/2 and Vinaxiz ~ 5000 calculated for (a) TEm cavity mode (following the dashed 
line) with higher modes giving lesser contributions, respectively; and (b) TMm cavity mode (following the dashed line), where 
TM112 < TM114 < TM113 at 500 ps, for example. Fig. 3 (a) inlay: Unitarity constraint, dm = 1 — X]n"°^(l^™"l^ ~ |/3mnP) in 
this case with m = 1 and cutoff imax = 51. For TMm the unitarity constraint is within 2 x lO^** (not shown here). 




FIG. 4: (Color online) Niii{t ~ IIT) for (a) TEm and (b) TMm modes at 77 = 1/2 with Vinaxiz = 5000. Open squares and 
filled circles denote single mode and multimopde {imax = 51) couplings, respectively, and arrows indicate optimum conditions. 



electric field component, E^, to the conduction electrons 
in the plasma. 

From an experimental viewpoint, it is desirable to de- 
crease the laser power for the suppression of temperature 
rises in the cavity. Shown in Figs. [l](c) and (d), respec- 
tively, are the variation of TEii„ and TMii„ modes, cal- 
culated assuming a very low laser power of O.Ol^J/pulse 
(Knax-^z = 1) at 77 = 1/2. Also iu this case, the plasma 
frequency of ~ 4 x 10^^ s~^ is still much larger than 
the frequency of TEii„ and TMii„ modes. As expected, 
the width and frequency variation of A:„c for TEii„ modes 
drastically decreases. Interestingly; however, such an ex- 
treme reduction of laser power does not decrease the fre- 
quency variation for TMii„ modes, although the widths 
are somewhat reduced. Such a large frequency variation 
leads to a large enhancement of photon creation and is 
intimately related to the TM plasma sheet boundary con- 
ditions. Indeed, as discussed later, the number of created 
photons for TMm docs not decrease even with a reduc- 



tion in laser power of order 1/5000 (0.01/iJ/pulsc), which 
is quite different from the TEm case. 

The time evolution of the number of photons for 
TEm/TMiii modes with SD position at the midpoint, 
1] ~ 1/2, is depicted in Figs. [21(a) and (b) respectively, up 
to a time lapse of IIT with a driving period of T = 105 
ps (frequency of 9.52 GHz) for TEm and T = 113 ps 
(frequency of 8.85 GHz) for TMm modes. Here, the sta- 
tionary (1, 1, 1) mode in the cavity (0.05 x 0.05 x 0.10 m'^) 
corresponds to a frequency of 4.5 GHz and thus, the driv- 
ing laser frequency is expected to be around 9.0 GHz 
(with period T = 111.11 ps). These optimum driving 
periods are correlated with shifted frequency fc„c cor- 
responding to half of the driving periods of 104.2 and 
116.1 ps for TEm and TMm modes, respectively, as can 
be derived from Figs. [2] (a) and (b). As we mentioned 
previously such a frequency shift comes from the inter- 
action between the cavity photons and conduction elec- 
trons. Other mode couplings (112,113, and 114) shown 
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FIG. 5: (Color online) Niii{t ~ IIT) for TEm and TMm 
modes dependent on SD location for a given optimum laser 
period T with X^naxiz = 5000. Filled triangles and circles 
denote a low laser power with Vmaxiz ~ 1- 

in Fig. [3] (a) are also slightly enhanced with time. Some 
intermode couplings satisfying Will +a;ii,i = i^^iii {j- hi- 
teger) enhance photon creation, while those far from the 
resonant condition damp the number of photons created. 
An important check of the numerical results is the unitar- 
ity constraint, Eq. (|2ip . which leads to an independent 
check of the convergence of the cut-off, imax — 51, in the 
mode sums. The inlay in Fig. [3] (a) shows the deviation 
from the unitarity condition, showing sufficient conver- 
gence (deviation: < 2 x 10"'''). Similar convergence is 
also confirmed for TMm modes (not shown here). 

The number of created photons iViii(t) for t ~ IIT de- 
pendent on driving laser period, T is shown in Fig. 2] for 
both TEui/TMiii modes at ?7 = 1/2 (midpoint). Inter- 
estingly, multimode coupling suppresses photon creation 
for TEiii modes, while it enhances the number of pho- 
tons for TMiii modes. Such an enhancement is caused 
by constructive (wm + = j^iii) mode couplings, 
as mentioned before and shown in Figs. 131(a) and (b). It 
should be emphasized that the optimum driving period, 
which is correlated with the frequency variation and in- 
termode couplings, is different from the resonant one cor- 
responding to the unperturbed field eigenfrequency (sta- 
tionary mode frequency) for all values of 77. The number 
of photons for the TEm mode takes a maximum value of 
4.3 at T 105 ps, which is almost 15 times that at 111.1 
ps {fl — 2cl>iii). In contrast, Nm for the TM case os- 
cillates periodically and is significantly smaller than the 
maximum of Nm for the TE case (at 77 = 1/2), because 
the higher modes such as fc2, ^3, •■• still have large fre- 
quency variation for TM modes, see Fig. [2] Importantly, 
however, the driving frequency of 17 = 2wiii (111.1 ps) 
still gives a large number of photon creation for TMm 
modes. Such a behavior is probably due to the differ- 
ent natures of TE modes (with component) and TM 
modes (with Ez component) that satisfy different jump 
conditions. The present results clearly show that para- 



metric photon creation takes place for TEm and TMm 
at 77 = 1/2 for a relatively wide range of driving periods. 
This is quite advantageous for the experimental detection 
scheme of DCE photons using a plasma sheet irradiated 
by a pulsed laser. 

The parametric photon creation rate depends on the 
location of the SD, rj = d/L^, as shown in Fig. [5] for 
TEm and TMm modes, respectively. Here, the driving 
period T is chosen to give the optimum value. The A^m 
value for TEm modes appears to increase drastically 
with increasing rj and reaches a maximum at 77 = 1/2 
(the midpoint), while a weak dependence on SD position 
is seen for TMm modes. We should note that Am at 
the midpoint for TEm is much larger (by an order of 
three) than that at 77 ^ 1 (near the cavity wall). This 
probably comes from the fact that the TEm mode has a 
belly (anti-node) at the midpoint (77 ~ 1/2). For 77 <C 1, 
A^iii(t) for TMm is much larger than that for TEm 
modes, as predicted thus far [y, [l^, [l^. As we previ- 
ously mentioned, a decrease in laser power does not re- 
duce the amplitude of kn for TMii„ modes, although it 
does for the TE case. Indeed, a decrease in laser power of 
the order of 1/5000 leads to a reduction of only half the 
number of photons for the TMm mode at the midpoint, 
as shown in Fig. [5] (filled triangle). Unfortunately; how- 
ever, a very low laser power reduces almost proportion- 
ately the number of Nm photons for TEm at 7; = 1/2 
and for TMm with 77 <C 1, as indicated in Fig. 5. As 
mentioned above, our numerical calculations without any 
perturbations show that the best location for the SD is 
the midpoint of the cavity for both TEm and TMm, in 
particular for TMm modes even with a low laser power. 



VI. CONCLUSION 

In conclusion our numerical results, without any per- 
turbations, indicate that the midpoint is the best location 
for a semiconductor diaphragm (SD) for both TEm and 
TMm modes, particularly for TM modes even with a 
low laser power of around 0.01 /iJ/pulse. It is empha- 
sized that such a low laser power suppresses heat losses 
and makes it possible to keep the cavity at a low tem- 
perature. The location of the midpoint for the SD is an 
advantageous configuration to suppress background sig- 
nals, because such a geometry enables us to separate the 
cavity into two parts, one of which is chosen as a photon- 
creation/detection chamber with the other used as a laser 
irradiation space. 

Another important point highlighted here are the fre- 
quency shifts mainly caused by the time-varying areal 
density of conduction electrons, ns{t), which are excited 
by the pulsed laser. This requires careful tuning of the 
driving laser frequency, although we do find a relatively 
wide range of driving period allowed to induce paramet- 
ric photon creation for both TEm and TMm modes. 
The present numerical analysis clearly shows that an SD 
irradiated by a pulsed laser behaves as a high-frequency 
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switching (transmission/reflection) actuator [3^, \^ in- 
stead of a mechanically vibrating cavity wall (mechani- 
cal oscillations with frequencies of GHz order seem tech- 
nically difficult to achieve). Also see [4l|, for other 
alternate experimental proposals. 

We have designed a new DCE detection system using 
highly-cxcitcd Rydberg atom (Rb or K) beams, acting 
as a single photon detector with a high sensitivity and 
a cavity made of niobium, which can be cooled down to 
~ 100 mK using a dilute refrigerator combined with a 
pulse tube and thus, we expect a high Q value better 
than 10^. (Note, Q-factors as high as 3 x 10^° have been 
reported in superconducting cavities [i^ HI] leading to a 
lifetime of order ~ 1 s, but the insertion of the SD will 
reduce this value.) The experimental set-up is now in 
progress. 



vector potential A [45[: 



K ^ 



m* 



-El 



(Al) 



Combining the continuity condition 4- V • K = with 
the Lorenz gauge condition then leads to 



dtAo 



(A2) 



where n is the bulk charge density and Aq is the scalar 
potential. As a result, the surface charge density, cr, is 
given by 



-Ao . 



(A3) 
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Appendix A: Plasma Sheet Model 

In this appendix we briefly discuss the plasma sheet 
model and derive the jump conditions for the TE and 
TM fields (*, $) given in Eqs. ([T]) and ©. The two- 
dimensional {{x,y) plane) current density K for conduc- 
tion electrons with relaxation time, t, under the presence 
of an external electric field E can be expressed using the 



It is convenient to now define two Hertz vectors Tie 
and n,„ as l2i,[2il: 



^0 = -V Ue 

e 



dU, 



dt 



V X n„ 



(A4) 



where e and fi are the bulk electric permittivity and mag- 
netic permeability of the SD, respectively. For electro- 
magnetic waves propagating along the z-axis of a rectan- 
gular cavity, the Hertz vector potentials can be reduced 
to the following ID expression 



n„, = * e. 



(A5) 



where Bz is a unit vector in the z-direction. If 5* (He) 
and $ (n.m) ai'c given the fields E and B are expressed 

as 



E 
B 



(A6) 



r 



a" 



where we have used the fact that V5_<I'(^') = 
^j^'($) = -k5^1'($). Note that $ and * satisfy the 

classical 3D wave equations [V^ - £[idl\^{^) 0. The 
boundary conditions for a stationary interface are, e.g., 
see fH: 



(D2-Di)-e, = a, 
X (H2-H1) = K, 



(B2 -Bi) -e^ = 
X (E2 - El) = . 



(A7) 



respectively. Substituting in our values of cr and K we 
find the following boundary interface conditions at z = d, 



9 ? 
disc $((i) = -/^ , 2 \ dz<^\z=d = -Mo ,, \ d:,<^\z=d 

disc B:, = k^disc *(d) = (A8) 

? ? 

G ill -■ G Tl 

disc dz^U^d - M— = MO— ^*(d) 

m* m* 

disc d,M.=d = 



where E and H are the electric and magnetic fields; and 
D and B are the electric and magnetic flux densities. 



The above boundary conditions come from the following 
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differential equations {^o — > 1) (A9) 



which are those presented in Eqs. © and ^ for TE and 



TM modes, respectively. 

1 e^n. 
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